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Today, it is possible to take a stroll through the 
high-dimensional state space of hydrodynamic 
turbulence and observe that turbulent trajecto- 
ries are guided by close passes to invariant solu- 
tions of the Navier-Stokes equations. Charting 
how close these passes are is a geometer's task, 
but in order to place them on a map, one first has 
to deal with families of solutions equivalent under 
the symmetries of a given flow. Evolution in time 
decomposes the state space into a 'spaghetti' of 
time trajectories. Continuous spatial symmetries 
foliate it like the layers of an onion. In this vi- 
sual tour of dynamics, we use a low-dimensional 
flow to illustrate how this tangle can be unraveled 
(symmetry reduction), and how to pick a single 
representative point for each trajectory (section 
it) and group orbit (slice it). Once the symme- 
try induced degeneracies are out of the way, one 
can identify and describe the prominent turbulent 
structures by a taxonomy of invariant building 
blocks (numerically exact solutions of the Navier- 
Stokes equations, finite sets of relative equilibria 
and infinite hierarchies of relative periodic orbits) 
and describe the dynamics in terms of near passes 
to their heteroclinic connections. 



I. INTRODUCTION 

Over the last decade, new insights into the dy- 
namics of moderate Re turbulent flows have been 
gained through visualizations of their oo-dimensional 
state spaces by means of dynamically invariant, represen- 
tation independent coordinate frames constructed from 
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physically prominent unstable coherent structures, here- 
after referred to as templates. Navigating and charting 
the geometry of these extremely high-dimensional state 
spaces necessitates a reexamination of two of the basic 
tools of the theory of dynamical systems: Poincare sec- 
tions and symmetry reduction. 

In quantum- mechanical calculations, one always starts 
out by making sure that the Hamiltonian has been 
brought to its symmetry-reduced block-diagonal, irre- 
ducible form; anything else would be sheer masochism. 
As the dynamical theory of turbulence is still in its in- 
fancy, symmetry reduction is not yet a common practice 
in processing turbulence data collected in experimental 
measurements and numerical simulations. Symmetry 
reduction of nonlinear flows is much trickier than the 
more familiar theory of irreducible representations for 
linear problems such as quantum mechanics, so most of 
our sketches illustrate the simplest case, the 1-parameter 
compact continuous group SO (2) symmetry. 

We show here how to bring the numerical or exper- 
imental data to a symmetry-reduced format before any 
further analysis of it takes place. Our tool of choice is the 
linear implementation of the method of slices. Here, we 
extend this local method to a global reduction of a turbu- 
lent flow by defining local 'charts', their borders, and the 
ridges that glue these linear tiles into an atlas that spans 
the ergodic state space region of interest. While 'charts' 
and 'atlases' are standard tools in geometry, the prescrip- 
tion for explicit construction of a symmetry-reduced state 
space presented here is, to our knowledge, new. We ex- 
plain the key geometrical ideas in simple but illustrative 
settings, eschewing the fluid dynamical and group theo- 
retical technicalities. 

Let us begin by defining a dynamical system comprised 
of a flow f l and the state space M. on which it acts. If a 
group G of continuous transformations acts on a continu- 
ous time flow, each state space point owns a set of tangent 
vectors (Fig. 1 (a)). Integrated in time, the velocity vec- 
tor v(x) traces out a trajectory f T (x) (Fig. 1 (b)). Apply- 
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FIG. 1. (a) In the presence of an iV-continuous parame- 
ter symmetry, each state space point x owns (N+l) tangent 
vectors: one v(x) along the time flow x(t), and the N group 
tangents £i(x), £2(2?), ••• ,tAr(x) along infinitesimal symme- 
try shifts, tangent to the N- dimensional group orbit M x . (b) 
Each point has a unique trajectory (blue) under time evolu- 
tion, (c) Each point also belongs to a group orbit (green) of 
symmetry-related points. For SO (2), this is topologically a 
circle. Any two points on a group orbit are physically equiv- 
alent, but may lie far from each other in state space, (d) 
Together, time-evolution and group actions trace out a wurst 
of physically equivalent solutions. 



ing the continuous transformations traces out a group or- 
bit M x = {gx | g G G} (Fig. 1 (c)). Together, time evo- 
lution and group actions trace out a complicated smooth 
manifold, hereafter affectionately referred to as a wurst 
(see Figs. 1 (d), 4(b) and 8), which we shall here teach 
you how to slice. 

A flow is said to have symmetry G if the form of 
evolution equations x = v(x) is left invariant, v(x) = 
g~ 1 v(gx) , by the set of transformations g G G. If a 
flow has symmetry, the simplest solutions are highly sym- 
metric invariant equilibria and relative equilibria studied 
in bifurcation-theory approaches to the onset of turbu- 
lence. Physicists love symmetries, 10 but nature often 
prefers solutions of no symmetry: while the flow equa- 
tions may be invariant under G, turbulent solutions are 
not. The highly symmetric solutions often lie far from the 
regions of state space explored by turbulence 11 and thus 
are of limited usefulness in understanding its dynamics. 
In contrast, the relative periodic orbits studied here are 
embedded in the turbulence, and capture its geometry 
and statistics. 

We can make headway in unraveling the tangle of 1- 
dimensional time trajectories with the notion of recur- 
rence. To quantify how close the state of the system at 
a given time is to a previously visited state, we need the 
notion of distance between two points in state space. The 
simplest (but far from the only, or the most natural) is 
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the Euclidean norm 

Hs _ x 'f = {x - x '\ x - x >) = J2(x - x't . (1) 

i=l 

For experimental data, a better norm, for example, might 
be a distance between digitized images. While in this 
paper we simply assume that a norm is given, its impor- 
tance cannot be overstated: the construction of invariant, 
PDE discretization independent state space coordinates, 
the symmetry reduction by minimization of the distance 
between group orbits undertaken in what follows, and the 
utility of the charts so constructed all depend on a well- 
chosen notion of distance in the high-dimensional state 
spaces we are charting here. 

Given a notion of distance, we can talk about a 'neigh- 
borhood', an open set of nearby states with qualitatively 
similar dynamics. Our main task in what follows will 
be to make this precise by defining a chart over a neigh- 
borhood and its borders. Given distances and neighbor- 
hoods, the next key notion is measure, or how likely a 
typical trajectory is to visit a given neighborhood. Af- 
ter some observations of a given turbulent flow, one can 
identify a set of templates, 6 points x'^\ j = 1,2, • • • in 
the state space representative of the most frequently re- 
visited features of the flow. 

Our goals here are two-fold: (i) In sect. II, we re- 
view the method of Poincare sections, with emphasis 
on two particular aspects that are applicable to high- 
dimensional flows: the construction of multiple local lin- 
ear charts and the determination of their borders, (ii) In 
sect. Ill, we discuss the effect of continuous symmetries 
on nonlinear flows, and in sect. IV we use the lessons 
learned from our discussion of Poincare sections to aid 
us in the reduction of continuous symmetries, and, thus, 
enable us to commence a systematic charting of the long- 
time dynamics of high-dimensional flows (sect. V). 

II. SECTION 

In the Poincare section method, one records the coor- 
dinates x n of the trajectory x(r) at the instants r n when 
it traverses a fixed oriented hypersurface V of codimen- 
sion 1. For the high-dimensional flows that we have in 
mind, the practical choice is a hyperplane, the only type 
of Poincare section (from now on, just a section) that we 
shall consider here. One can choose a section such that it 
contains a template of interest. Properly oriented, such 
a section can capture important features of the flow in 
the neighborhood of the section- fixing template. 

But how far does this neighborhood extend? The an- 
swer is that the section captures neighboring trajectories 
as long as it cuts them transversally; it fails the moment 
the velocity field at a point x* fails to pierce the section. 
At these locations, the velocity either vanishes (equilib- 
rium) or is orthogonal to the section normal n, 

n-v(x*) = 0, x*eS. (2) 
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For a smooth flow in d dimensions such points form a 
smooth (d — 2)-dimensional section border 5 c?, which 
encloses the open neighborhood of the template charac- 
terized by qualitatively similar flow. We shall refer to this 
region of the section as a chart of the template neighbor- 
hood (see Fig. 2). Beyond the border, the flow pierces 
the section in the 'wrong' direction and the dynamics are 
qualitatively different. 

As an example consider the Rossler system , 

x — —y — z 

y = x + ay (3) 
z = b + z(x — c) , 

where a = b = 0.2 and c = 5.7. This flow has two 
prominent invariant states, the 'inner' and the 'outer' 
unstable equilibria x'^ and x'^ (see Fig. 2 (a)) , which 
we choose as templates for our sections. 

We orient the sections so the plane V- contains x'^ 
and its 1-dimensional stable eigenvector (Fig. 2 (b)), and 
the other section T+ contains x'^ and its 1-dimensional 
unstable eigenvector (Fig. 2(c)), thus capturing the lo- 
cal spiral-in, spiral-out dynamics. The remaining free- 
dom to rotate each section can be used to orient them 
in such a way that the ridge (the intersection of the two 
sections) lies approximately between the two templates 
(Fig. 2 (d)). Choosing sections is a dark art: in the exam- 
ple at hand the dynamics of interest is captured by the 
two charts - if that were not the case, one would have had 
to interpolate, by inserting a third chart between them. 

For Rossler flow, the border condition (2) yields a 
quadratic condition in 3 dimensions, so the section bor- 
ders drawn in Fig. 2 (b) and Fig. 2 (c) are conic sections. 
The two charts meet at a ridge, and together do a pretty 
good job as the 2-chart atlas of the interesting Rossler 
dynamics. Due to the extreme contraction rate of the at- 
tractor, its intersection with the section in Fig. 2 (b) is for 
all practical purposes 1-dimensional, and the associated 
return map yields all periodic orbits of the 3-dimensional 
flow. 13 

In 3 dimensions everything — sections, ridges, sec- 
tion borders — can be drawn and the chart fits on a 
2-dimensional sheet of papyrus. But what about for 
hydrodynamic flows where the dimensionality d of the 
state space is very large? The point of the cartographical 
enterprize undertaken here is that while it is impossible 
to visualize the (d— 2)-dimensional section border of the 
(d—1) -dimensional slab that is now our chart, a point 
is a point and a line is a line in a projection from any 
number of dimensions, so a trajectory crossing of either a 
section or a section border can be easily determined and 
visualized in any dimension. 

To summarize: Evolution in time decomposes the state 
space into a spaghetti of 1-dimensional trajectories x(r), 
each fixed by picking a single point x(0) on it. A well 
chosen set of section charts of codimension 1 allows us to 
'quotient' the continuous time parameter r, and reveal 
the dynamically important transverse structure of the 
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FIG. 2. 2-chart atlas for Rossler flow, (a) The inner equilib- 
rium x'^ is a (spiral-out) saddle-focus with a 2-dimensional 
unstable manifold and a 1-dimensional stable manifold. The 
outer equilibrium x'^ is a (spiral-in) saddle- focus, with a 2- 
dimensional stable manifold (basin boundary for initial con- 
ditions that either fall into the chaotic attractor, or escape to 
infinity) and a 1-dimensional unstable manifold, (b) Chart 
V- of the x'^ neighborhood carved out of a Poincare sec- 
tion plane through the inner equilibrium x'^ and its stable 
eigenvector, with section border drawn as the solid red line. 
Note the ridge (dashed blue line): the chart stops at the ridge, 
(c) Chart V+ (here viewed from below) is bounded by section 
border (solid red line) of a section through the outer equi- 
librium x'^ and its unstable eigenvector. The chart stops 
at the ridge (dashed blue line), and it does not intersect the 
strange attractor. (d) A two-chart atlas of Rossler flow, with 
charts V- and V+ oriented and combined so that the ridge 
(intersection of the two sections, indicated by the dashed blue 
line in the three figures) lies approximately between the tem- 
plates. Section hyperplanes beyond this ridge do not belong 
to the atlas. 



flow's stable/unstable manifolds. For unstable trajecto- 
ries one needs, in addition, a notion of recurrence to the 
section. The set of points {x n } = {x(r n )}, separated by 
short time flights in between sections, captures the trans- 
verse dynamics without losing any information about the 
chaotic flow. We can thus chart interesting regions of 
state space by picking a sufficient number of templates 
and using them to construct charts of their neighbor- 
hoods, each bounded by section borders and ridges. 

We close this section with a remark on what sections 
are not: A Poincare section is not a projection onto 
a lower-dimensional space (in sense that a photograph 
is a 2-dimensional projection of a 3-dimensional space). 
Rather, it is a local change of coordinates to a direction 
along the flow v(x), and the remaining coordinates trans- 
verse to it. No information about the flow is lost; the full 
space trajectory x(r) can always be reconstructed by in- 
tegration from its point x in the section. 
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dynamical system with a continuous (but no discrete) 
symmetry, equivariant under SO(2) rotations by 



FIG. 3. A symmetry relates physically equivalent states; a 
pipe flow solution translated or rotated is also a solution, (a) 
An instantaneous state of the fluid is indicated by a 'swirl' 
- here the reader has to imagine a particular instantaneous 
velocity field across the entire pipe. The same state may be 
rigidly (b) translated by downstream shift £ (fluid states are 
SO (2) z equivariant in a stream-wise periodic pipe), (c) trans- 
lated by t and rotated azimuthally by (the two states are 
SO (2)e x SO (2) z equivariant), and (d) reflected and rotated 
azimuthally by <j) (the two states are 0(2)q equivariant). Some 
symmetry-related states may also be connected by time evo- 
lution. A relative equilibrium is a solution of the equations 
of motion that retains its shape while rotating and traveling 
downstream with constant phase velocity c. A relative peri- 
odic orbit M p is a time dependent, shape-changing state of 
the fluid that after a period T p reemerges as (b), (c), or (d), 
the initial state translated by £ p , rotated by <j) p and possibly 
also azimuthally reflected. 



III. DANCERS AND DRIFTERS 

What is a symmetry? A visualization of the fluid dy- 
namics of a pipe flow, Fig. 3, affords an intuitive illustra- 
tion. Solutions of pipe flow remain physically the same 
under azimuthal rotations and stream-wise translations 
(which become SO (2) rotations in numerical stream- wise 
periodic pipes), but rotated and shifted solutions may 
correspond to distant points in state space. 

Each SO (2) group orbit is topologically a circle, but 
it traces out a complicated state space curve composed 
of many Fourier modes that are nonlinearly coupled and 
thus of comparable magnitude. Together, the two SO (2) 
symmetries of numerical pipe flow sweep out contorted 
and hard to visualize T 2 tori (see ref. ), so we shall 
illustrate the key ideas by a much simpler example, the 
S0(2)-equivariant Gibbon and McGuinness ' complex 
Lorenz equations of geophysics and laser physics, 



x\ = — axi + cryi , x 2 = —crx 2 
?/i = (pi - z)x 1 - p 2 x 2 -yi~ ey 2 
2/2 = P2X1 + (pi ~ z)x 2 + ey 1 - y 2 
z = -bz + xiyi + x 2 y 2 • 



cry 2 



(4) 



Here all coordinates and parameters are real. In our cal- 
culations the parameters are set to those used by Simi- 
nos' p x = 28, p 2 = 0, b = 8/3, a = 10, e = 1/10. The 
complex Lorenz equations are an example of a simple 
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The group is 1-dimensional and compact, its elements 
parameterized by (j) mod 2tt. For historical background, 
Poincare return maps, symbolic dynamics and in-depth 
investigation of the model, see refs. 8 ' 17 . 

The strange attractor of the complex Lorenz flow , in 
its present state, is a complete mess (Fig. 4(a)). Solu- 
tions tend to drift along continuous symmetry directions, 
with the physically important shape-changing dynamics 
hidden from view. 

The ultimate drifter, the signature invariant solution 
that signals the presence of a continuous symmetry is 
a relative equilibrium (traveling wave, rotational wave, 
etc.), a trajectory whose velocity field lies within the 
group tangent space, v(x) = c-t(x) , and whose time evo- 
lution is thus confined to the group orbit (see Fig. 4 (b)); 
think of an unchanging body carried by a stream. 

A relative periodic orbit behaves more like a dancer. 
M p is a trajectory that recurs exactly 



x(t) = g p x(t + T p ) , x(r) e M p , 



(5) 



after a fixed relative period T p , but shifted by a fixed 
group action g p that maps the endpoint x(T p ) back into 
the initial point cycle point x(0); think of a dancer mov- 
ing across the stage through a set of motions and then 
striking her initial pose, or study the pipe flow sketches 
in Fig. 3. 

Because the SO (2) transformations act on the complex 
Lorenz flow only through the simplest m = 1 Fourier 
mode, here all group orbits are circles and appear ellip- 
tical in d = 5 — >• 3 dimensions projections. Nevertheless, 
even the wurst traced out by one of the simplest relative 
periodic orbits 01 (shown in Fig. 4 (b)) is not so easy to 
get one's head around: you are looking at a 3-dimensional 
projection of a torus embedded in 5 dimensions. 

To summarize: continuous symmetries in the dynam- 
ics foliate the state space into an onion, where each layer 
is a group orbit (Fig. 5(a)). How are we to sort out 
this mess? All the points on a group orbit are physi- 
cally equivalent, so we are free to replace a given flow 
x(t) by any other x(r), such that x(r) — g(r) x{r) by 
a moving frame 19 ' 21 transformation g(r). As long as no 
symmetry reduction procedure is prescribed, g(r) is free: 
it can be any, in general time dependent, group transfor- 
mation. For example, to film our dancer, we can mount 
the camera on a cart moving alongside her. So, in the 
presence of continuous symmetries, there are two kinds of 
motion: those of a dancer, continuously changing shapes, 
and those of a drifter, merely shuffling along the shape 
invariant directions. We will presently banish the drifters 
and just enjoy the dance. 
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FIG. 4. Complex Lorenz flow, d = 5 —> 3 dimensional 
{x±,X2,z} projections: (a) The strange attractor. (b) The 
initial relative equilibrium TW\ point is shown by the red 
dot, and its group orbit / trajectory by the dashed red line. 
One period of the 01 relative periodic orbit is shown by the 
solid blue line. The group orbit of its (arbitrary) starting 
point is shown by the dashed blue line: after one period the 
trajectory has returned to the group orbit but with a different 
phase. The wurst, i.e., the group orbit of the 01 trajectory 
(dark blue) is shown by the cyan surface. Following 01 for 15 
more periods (faint dotted lines) starts filling out this torus; 
in that time, the slowly drifting relative equilibrium TW\ has 
advanced to the next red dot (red line). Symmetry-reduced 
complex Lorenz flow, d = 4 — » 3 dimensional {x2,2/2, z} pro- 
jections: (c) Strange attractor from frame (a) reduced to a 
single slice hyperplane, using TWi as the template. 01 is 
now a periodic orbit, shown by the solid black line. The 
dynamics exhibits singular jumps (shown in red) due to for- 
bidden crossings of the chart border. In contrast to the 1- 
dimensional section borders of Fig. 2, here the chart borders 
are 3-dimensional and hard to visualize, (d) The 2-chart at- 
las (see the sketch of Fig. 11) of the same strange attrac- 
tor encounters no chart borders and exhibits no singularities. 
The trajectory changes colors from red to blue as it crosses 
between the slice hyperplanes of x'^ and x'^ 2 \ The ridge 
(shown in brown) acts as a Poincare section V with red or 
blue ridge points x* marking the direction of the crossing. 
The charts are 4-dimensional, the ridge 3-dimensional, so the 
colored blocks and planes are only cartoon drawings of their 
projections onto the 2-dimensional figure. 



IV. CHART 

Suppose you are computing a set of numerically exact 
invariant solutions of the Navier- Stokes equations. Do 
you want to compute the same solution over and over 
again, once for every point on the group orbit? No, you 
would like to compute it only once. The strategy for pick- 
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FIG. 5. (a) The iV-dimensional group orbit M x (o) of state 
space point x(0) and the group orbit M X ( T ) reached by the 
trajectory x(r) a time r later, (b) Two physically equiva- 
lent trajectories x(r) and x(r) are related, in general, by an 
arbitrary, time dependent moving frame transformation g(r), 
such that x(t) = g(r) x(r). (c) A symmetry reduction scheme 
M — »• M is a rule that prescribes g(r) and thus replaces a 
group orbit M x C M through x by a single point x £ M. 
(d) In this paper, g(r) is fixed variationally by the extremal 
condition (7) for the point x on the group orbit M x that is 
nearest to the template x! . 

ing out that one representative solution is called symme- 
try reduction. Its goal is to replace each group orbit by 
a unique point in a lower-dimensional symmetry-reduced 
state space M c M/G, as sketched in Fig. 5 (c). 

What is a smart way to go about it? Intuition gained 
from pipe flow (see Fig. 3) will again prove helpful. A 
turbulent flow exhibits a myriad of unstable structures, 
each traveling down the pipe with its own phase veloc- 
ity. The method of slices that we now describe tells 
you how to pull each solution back into a fixed frame 
called a slice and compare it to your repertoire of pre- 
computed solutions, or the templates {x'^}, using the 
poor geometer's version of a geodesic, the principle of 
the closest distance to each. What follows is similar to 
the construction of sections of sect. II; due to the linear 
action of the symmetry group, slicing is easier than sec- 
tioning, but wholly unfamiliar. This is why we reviewed 
the Poincare sections first. We now offer a pictorial tour 
of this (save for one bold incursion ) hitherto uncharted 
territory. 

First, pick a template x 1 and use the freedom to shift 
and rotate it (Fig. 5 (b)) until it overlies, as well as pos- 
sible, the state x, by minimizing the distance 

\\x-g{4>)x'\\. (6) 

Now, replace the entire group orbit of x by the closest 
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match to the template pattern, given by x = g~ x x. From 
here on, we will use the hat on x to indicate the unique 
point on the group orbit of x that is closest to the tem- 
plate x' . The symmetry-reduced state space M. is com- 
prised of such closest matches, a point for each full state 
space group orbit. 

The minimal distance satisfies the extremum condition 
(Fig. 5 (d)) 

- g(<j>) x'f =2(x- x'\t') = , t' = Tx', 

ocp 

where the [d x d] matrix T is the generator of infinites- 
imal symmetry transformations. = is a 
constant. To streamline the exposition, we shall assume 
here that the symmetry group is SO(n). In that case T is 
antisymmetric, so the group tangent vector t' evaluated 
at x! is normal to x! and the term (x'\Tx') vanishes. 
Therefore x , the point on the group orbit of x that lands 
in the slice satisfies the slice condition 

(x\t')=0. (7) 

As x{t) varies in time, the template x' tracks the 
motion using the slice condition (7) to minimize 
\\x(r) — g(4>(r))x f \\^ and the full-space trajectory x(r) is 
rotated into the reduced state space trajectory x(r) by 
appropriate time varying moving frame angles {(/>(r)}, as 
depicted in Fig. 6(a). M. is thus a [d — N) -dimensional 
hyperplane normal to the N group tangents evaluated at 
the x' as sketched in Fig. 6 in a highly idealized manner: 
A group orbit is an Af-dimensional manifold and, even 
for SO (2), is usually only topologically a circle and can 
intersect a hyperplane any number of times (see Figs. 7 
and 8). 

One can write the equations for the flow in the reduced 
state space x = v{x ) (for details see, for example, ref. 13 ) 
as 

v(x) = v(x) — <j)(x)t(x) (8) 

m = (v(x)\t')/(tm'), (q) 

which confines the motion to the slice hyperplane. Thus, 
the dynamical system /*} with continuous symme- 
try G is replaced by the reduced state space dynamics 
/*}: The velocity in the full state space v is the 
sum of v, the velocity component in the slice hyperplane, 
and 0t, the velocity component along the group tangent 
space. The integral of the reconstruction equation for <\> 
keeps track of the group shift in the full state space. 

The template x 1 should be a generic state space point 
in the sense that its group orbit has the full N dimen- 
sions of the group G. The set of the group orbit points 
closest to the template x 1 forms a neighborhood of x' 
in which each group orbit intersects the hyperplane only 
once. A slice hyperplane qualitatively captures neighbor- 
ing group orbits until, for a point x* not so close to the 
template, the group tangent vector t(x*) lies in the slice 
hyperplane. The group orbits for such points are grazed 
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FIG. 6. The method of slices, a state space visualiza- 
tion: (a) A chart M C M/G lies in the (d— A/") -dimensional 
slice hyperplane (7) normal to ti...t' N , which span the N- 
dimensional space tangent to the group orbit gx (dotted 
line) evaluated at the template point x . The hyperplane in- 
tersects all full state space group orbits (green dashes). The 
full state space trajectory x(r) (blue) and the reduced state 
space trajectory x(r) (green) are equivalent up to a 'moving 
frame' rotation x(r) — g(r)x(r), where g(r) is a shorthand 
for g((j>(r)). (b) In the full state space, a relative periodic 
orbit x(0) —> x(r) — ► x(T p ) returns to the group orbit of x(0) 
after a time T p , such that x(0) = g p x(T p ). A generic relative 
periodic orbit quasi-periodically fills out what is topologically 
a torus (Fig. 4(b)). In the slice, the symmetry-reduced tra- 
jectory is periodic, x(0) = x(T p ). 

tangentially rather than sliced transversally, much like 
what happens at the section border (2) for evolution in 
time. This is also a linear condition and defines the chart 
border <S, ' a (d — N— l)-dimensional manifold, which 
contains all the points x* whose group tangents lie in the 
slice hyperplane, i.e., 

(£*!£') = and |0 = 0. (10) 

S also contains all points for which t(x*) = 0. While for 
the Poincare sections (2) the analogous points were equi- 
libria (captured only if the section cut through them), 
for slice hyperplanes points with vanishing group actions 
belong to invariant subspaces, and, by its definition, ev- 
ery chart border automatically includes all invariant sub- 
spaces. 

For the complex Lorenz equations (4), the invariant 
subspace is the 1-dimensional z-axis, with trivial dynam- 



High-dimensional cartography 



7 




FIG. 7. The chart border is the (d—N— l)-dimensional hy- 
perplane that contains all the points x* whose group tangents 
t(x*) lie in the slice hyperplane or vanish and are thus normal 
to t' . Beyond this boundary, the group orbits pierce the slice 
hyperplane in the wrong direction, so only the half- hyperplane 
that contains the template belongs to the slice. The chart bor- 
der is not easy to visualize; For the lack of dimensions, here 
it is drawn as a 'line', the z axis in this 3-dimensional sketch, 
(a) If the equivariant coordinates transform only under the 
m = 1 representation of SO(2), every group orbit is a circle, 
and crosses any slice hyperplane exactly twice. However, if 
there are coordinates that transform as higher m, the group 
orbit can pierce the hyperplane up to 2m times, and the chart 
border lies closer to the template: For example, (b) a group 
orbit for a combination of m = 1 and m = 2 equivariant co- 
ordinates resembles the seam of a baseball, and can cross the 
slice hyperplane 4 times, out of which only the point closest 
to the template is in the slice. 




FIG. 8. Wurst, sliced. Every slice hyperplane cuts every 
group orbit at least twice (see Fig. 6), once at the orbit's 
closest passage to the template, and another time at the most 
distant passage, also satisfying the slice condition (7). An 
SO (2) relative periodic orbit M p is topologically a torus, so 
the two cuts are the two periodic orbit images of the same 
relative periodic orbit, the good close one x p (blue), and the 
bad distant one (red), on the other side of chart border, and 
thus not in the slice. 



ics, z = —bz, but in general invariant subspaces are 
high-dimensional and have their own dynamics. Physi- 
cists, for example general relativists, often work in in- 
variant subspaces, as this is easier than solving the full 
problem. Such approaches yield highly symmetric so- 
lutions, 23 ' 24 whose dynamics may be quite different from 
those that guide turbulence in the full state space (for a 
striking example, see ref. 11 ). 



There is yet another, much kinder type of a border: 
a ridge. Our initial chart Ai^ is a (d— TV) -dimensional 
hyperplane. If we pick another template point x'^ 2 \ it 
comes along with its own slice hyperplane A4^ 2 \ Any 
pair of (d— 7V)-dimensional local slice hyperplanes inter- 
sects in a ridge, a (d—N — l)-dimensional hyperplane V 
of points shared by a pair of charts and thus satisfying 
the slice condition (7) for both, 

{x*\1/W) = and (x*\1/W) = 0. (11) 

The ridge forms a Poincare section that serves as 

a toll bridge, crossed by any direct transit from a chart 
to the adjacent chart In Fig. 10(a) a ridge 

is visualized as a 'line', and in Fig. 11 as a 'plane' of 
intersection of two volumes. We shall refer to the neigh- 
borhood of a template x'^ bounded by its chart border 
and the ridges to other such linear neighborhoods as a 
chart M^ j) C M/G, and to (10) and (11) as the border 
conditions. 



V. CHARTING THE SLICE 

Let us summarize the voyage so far: we are charting 
a curved manifold, and it would be nice to use tools 
of differential geometry, but this seems not possible in 
the high-dimensional state space of hydrodynamics tur- 
bulence. The only feasible way to chart this space is to 
(1) quotient all continuous symmetries, and (2) tile the 
reduced state space with flat (d — 7V)-dimensional tiles, 
or charts. We do this step by step, starting with a set 
of templates and using them to construct charts of each 
neighborhood, and then building up an atlas of the slice, 
chart by chart, which captures all of the reduced dynam- 
ics of interest (but not all possible dynamics). Here are 
the steps along the way: 

Template: Pick a template x ' such that G acts on it 
regularly with a group orbit of dimension N . 

Slice hyperplane: The (d— 7V)-dimensional hyperplane 
satisfying (x\t f a ) = , normal to group transforma- 
tion directions at the template x' . 

Moving frame: For any x, the slice condition (x\t f ) = 
on x = g(4>)x determines the moving frame, i.e., 
the group action g(<j)) that brings x into the slice 
hyperplane. 

Chart border: The set of points x* on a slice hyper- 
plane whose group orbits graze the hyperplane tan- 
gentially, such that (x*\t f ) = (t(x*)\t') = 0. 

Flow invariant subspace: If a subset or all of the 
group tangents of a chart border point x* vanish, 
ta( x *) = 0, its time trajectory remains within a 
flow- invariant subspace for all times. 
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Ridge: A hyperplane of points x* G V^ 21 ^ formed by the 
intersection of a pair of slice hyperplanes Ai^ and 

Chart: The neighborhood of a template x'^\ bounded 
by the chart border and the ridges to other linear 
neighborhoods, comprises a chart jCi^ c M/G. 
The borders ensure that there is no more than one 
oriented group orbit traversal per chart; a group 
orbit either pierces one chart, or no charts at all. 

Atlas: A set of (d — A^)-dimensional contiguous charts 

Slice: Let G act on a d-dimensional manifold M, with 
group orbits of dimension N or less. A slice is a 
(d — 7V)-dimensional submanifold M such that all 
group orbits that intersect M do so transversally 
and only once. 

In the literature, 25-27 'slice' refers to any co- 
dimension N manifold that slices transversally a group 
orbit. Here, we define an atlas over a slice constructively 
but more narrowly, as a contiguous set of flat charts, 
with every group orbit accounted for by the atlas sliced 
only once, and belonging to a single chart. A slice is not 
global, it slices only the group orbits in an open neigh- 
borhood of the state space region of interest. 

The physical task, for a given dynamical flow, is to pick 
a set of qualitatively distinct templates (for a turbulent 
pipe flow there might be one typical of 2-roll states, one 
for 4-roll states, and so on), which together provide a 
good atlas for the region of M/G explored by chaotic 
trajectories. 

The rest is geometry of hyperplanes and has nothing 
to do with dynamics. Group orbits M x u) through x^\ 
group tangents t(x'^), and the associated charts M^ 
are purely group-theoretic concepts. The slice, chart bor- 
der and ridge conditions (7), (10) and (11) are all linear 
conditions which depend on the ray defined by the tem- 
plate x', not its magnitude. Checking whether the chart 
border is on the far side of the ridge between two slice 
hyperplanes is a linear computation; for a symmetry- 
reduced trajectory moving in M^ chart one only has 
to keep checking the sign of 

(x(r)\t'W). (12) 

Once the sign changes, the ridge has been crossed, and 
from then on the trajectory should be reduced to the 
M^ chart. For three or more charts you will have to 
align the ridge of the current chart with a previously- 
used chart. You'll cross that ridge when you come to it 
(a hint: the manifold is curved, so there will be a finite 
jump in phase). 

How the charts are put together is best told as a 
graphic tale, in the 5 frames of Figs. 9, 10 and 11, and 
then illustrated by contrasting the mess of the complex 
Lorenz equations strange attractor Fig. 4 (a) to the ele- 
gance of its 2-chart atlas, Fig. 4 (d). 
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(a) 



(b) 




FIG. 9. A 2-chart atlas. Sketch (a) depicts two templates 
x'^ 2 \ each with its group orbit. Start with the template 
x^. All group orbits traverse its (d — l)-dimensional slice 
hyperplane, including the group orbit of the second template 
x'^ 2 \ (b) Replace the second template by its closest group 
orbit point x'^ 2 \ i.e., the point in chart M^\ This is allowed 
as long as x'^ is closer than the M^ chart border (red 
region), otherwise an interpolating, closer template needs to 
be introduced. 

It is worthwhile to note that the only object that enters 
the slice hyperplane, border and ridge conditions is the 
ray defined by the unit vector i = t'/\\t'\\. This gives 
much freedom in picking templates. In particular, the 
two rays 

= (0.263,-0.692,0.624,-0.251,0) 
f /(2) = (0.153, -0.610, 0.747, -0.213, 0) (13) 

used to construct the complex Lorenz equations 2-chart 
atlas of Fig. 4 (d) were found by numerical experimenta- 
tion. 

With the atlas in hand, the dynamics is fully charted: 
as explained in refs. 8 ' 13 , Poincare return maps then yield 
all admissible relative periodic orbits. 

Three concluding remarks on what slices are not: 
(1) Symmetry reduction is not a dimensional-reduction 
scheme, a projection onto fewer coordinates, or flow mod- 
eling by fewer degrees of freedom: It is a local change of 
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: '(2) Q 



(a) 




(b) 



FIG. 10. A 2- chart atlas, (a) Now that the group orbits have 
been reduced to points, erase them and consider the two slice 
hyperplanes through the two templates. As these two tem- 
plates are the closest points viewed from either group orbit, 
they lie in both slice hyperplanes. However, the two tangent 
vectors t'^ and t ,(y2 ^ have different orientations, so they de- 
fine two distinct charts jCi^ and which intersect in the 
ridge, a hyperplane of dimension (d—2) (here drawn as a 'line', 
and in Fig. 11 as intersection of two 'volumes') shared by the 
template pair that satisfies both slice conditions (11). The 
chart for the neighborhood of each template (a page of the 
atlas in part (b)) extends only as far as this ridge. If the 
templates are sufficiently close, the chart border of each slice 
hyperplane (red region) is beyond this ridge, and not encoun- 
tered by the symmetry-reduced trajectory x(r). The reduced 
trajectory is continuous in the slice comprised of such charts 
- it switches the chart whenever it crosses a ridge, (b) The 
slice (unique point for each group orbit) is now covered by an 
atlas consisting of (d— l)-dimensional charts M^ 2 \ 



coordinates with one (or N) coordinate (s) pointing along 
the continuous symmetry directions. No information is 
lost by symmetry 'reduction', one can go freely between 
solutions in the full and reduced state spaces by integrat- 
ing the associated reconstruction equations (9). 

(2) If the flow is also invariant under discrete symme- 
tries, these should be reduced by methods described, for 
example, in ChaosBook.org. 

(3) An atlas is not needed for Newton determination 
of a single invariant solution, or a study of its bifurca- 
tions. Any local section and slice plus time and shift 
constraints does the job. It is possible to compute 




FIG. 11. Here the two charts of Fig. 10 (a) are drawn as two 
(d— l)-dimensional slabs. The ridge, their (d— 2)-dimensional 
intersection, can then be drawn as the shaded plane. This 
hyperplane cuts across the symmetry-reduced trajectory x(r) 
and thus serves as a Poincare section V^ 2V) that captures all 
transits from the neighborhood of template to the neigh- 
borhood of template x^ 2 \ Poincare section transits are ori- 
ented, so x\ and X2 are in the section, but the third point is 
not. 



60,000 relative periodic orbits this way. 22 Once we have 
more than one invariant solution, the question is: how is 
this totality of solutions interrelated? For that, a good 
atlas is a necessity. 



VI. BRIDGES TO NOWHERE 

Everybody encounters a symmetry sooner or later, so 
the literature on symmetry reduction is vast (for a histor- 
ical overview, see remarks in ChaosBook.org and ref. 8 ). 
Before asking, "Why the method of slices and not [...]?" 
a brief tour of the more familiar symmetry reduction 
schemes is called for. They all have one thing in com- 
mon: they will not work for high-dimensional nonlinear 
systems. 

To start with, mastery of quantum- mechanics or bi- 
furcation theory 28 ' 32 symmetry reduction to linear irre- 
ducible representations is only partially illuminating; lin- 
ear theory works quite well for linear unitary operators 
or close to a bifurcation, but, as we tried to show in 
this pictorial tour, the way symmetries act on nonlinear 
systems is much subtler. For flows with strongly nonlin- 
early coupled modes, both time trajectories and group 
orbits are complicated, so choices of sections and slices 
require insight into the geometry of the particular flow, 
there exists no general theory of linear transformations 
into symmetry irreducible coordinates that would do the 
job. 

There are purely group-theoretical approaches, with 
no dynamics to inform them, inspired by the observation 
that while coordinates xi are equivariant, the squared 
length r 2 = ^2 x\ is invariant under 0(n) transforma- 
tions. 
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For SO (2), an obvious idea is to go to polar coordi- 
nates. The simplest nonlinear examples already run 
into Tj — >• type of singularities, and it is not alto- 
gether clear how one would rewrite the Navier- Stokes 
equations in such a format, or integrate them numeri- 
cally. A more sophisticated approach is to rewrite the 
dynamics in terms of invariant polynomial bases, de- 
scribed lucidly in ref. 34 , with the equivariant state space 
coordinates (xi, #2, #3, xj) replaced by an invariant 
polynomial basis (wi,W2,W3,...,w m ). As the dimension of 
the problem increases, the number of these polynomials 
grows quickly, as does the number of syzygies, the non- 
linear relations amongst them. There is no guiding prin- 
ciple for picking a set of such polynomials, and no practi- 
cal way to implement the scheme for high-dimensional 
flows: how and why would one replace the large number 
of equivariant state space coordinates of hydrodynamic 
turbulence with a vast number of invariant polynomials? 

Others approaches are informed by dynamics, foremost 
among them being the method of co-moving frames. Vi- 
sualizing a single 'relative' trajectory in its co-moving 
frame, i.e., moving with that solution's mean phase ve- 
locity, is useful if one is concerned with that individual 
solution and the tiny relative periodic orbits (modulated- 
amplitude waves) that bifurcate off it. ' A co- moving 
frame is useless, however, if we are concerned with study- 
ing collections of these trajectories, as each solution trav- 
els with its own mean phase velocity c p = <fi p /T p , and 
there is no single co-moving frame that can simultane- 
ously reduce all traveling solutions. The slice that we 
construct here is not 'co- moving', but emphatically sta- 
tionary. 

There exists a beautiful theory of symplectic symme- 
try reduction for the mechanics of three-dimensional rigid 
bodies, 36 ' 37 or using Lie symmetry reduction to derive 
Eulerian velocity fields from Lagrangian trajectories. 38 
These approaches do not appear to be applicable to prob- 
lems considered here, and anyway, the goal is different. 
Rather than to reduce a particular set of equations, we 
seek to formulate a computationally straightforward and 
general method of reducing any continuous symmetry, 
for any high-dimensional chaotic/turbulent flow. One 
should also note that 'symmetry reduction' in general 
relativity 23 and Lie theory often implies restricting one's 
solution space to a subspace of higher symmetry; here we 
always work in the full state space. 

There is, however, one intriguing, compelling and phys- 
ically informed contender. In mechanics and field theory 
it is natural to separate the flow locally into group dy- 
namics and a transverse, 'horizontal' flow, 37 ' 39 by the 
'method of connections', 40 illustrated in Fig. 12. The 
method of connections, however, does not reduce the dy- 
namics to a lower-dimensional reduced state space M/G. 
In contrast to the method of co-moving frames, where one 
defines a mean phase velocity of a relative periodic orbit, 
the method of connections is inherently local. The two 
methods coincide for relative equilibria. 

The meaning of the 'method of connections' in clas- 



10 




FIG. 12. (a) By equivariance v(x) can be replaced by v±(x), 
the velocity normal to the group tangent directions at state 
space point x. (b) The method of connections replaces v(x) 
at every instant x — x(r) by v±(x), so in x(r)'s covariant 
frame there is no motion along the group tangent directions. 

sical dynamics is clearest in the work of Shapere and 
Wilczek: 18 ' 41 one can observe a swimmer (or our dancer) 
from a fixed slice frame, or bring her back to observe only 
the shape-changing dynamics, no drifting. Left to her- 
self, she will reemerge in the same pose someplace else: 
that shift is called a 'geometrical phase', which -while ac- 
cruing it is the whole point of swimming- has not played 
any role in our discussion of symmetry reduction. Con- 
versely, most gauge choices in quantum field theory are 
covariant, and while that suffices to regularize path inte- 
grals, the method of slices says that this is no symmetry 
reduction at all, and it yields no insight into the geometry 
of nonlinear flows. 

Symmetry reduction in dynamics (including classical 
field theories such as the Navier-Stokes equations) closely 
parallels the reduction of gauge symmetry in quantum 
field theories. There, the freedom of choosing moving 
frames shown in Fig. 5 is called 'gauge freedom' and a 
particular prescription for choosing a representative from 
each gauge orbit is called 'gauge fixing'. Just like the slice 
hyperplanes of Fig. 7 may intersect a group orbit many 
times, a gauge fixing submanifold may not intersect a 
gauge orbit, or it may intersect it more than once ('Gri- 
bov ambiguity'). ' In this context a chart is called a 
'Gribov' or 'fundamental modular' region and its border 
is called a 'Gribov horizon' (a convex manifold in the 
space of gauge fields). The Gribov region is compact and 
bounded by the Gribov horizon. Within a Gribov region 
the 'Faddeev-Popov operator' (analogue of the group or- 
bit tangent vector) is strictly positive, while on the Gri- 
bov horizon it has at least one vanishing eigenvalue. 

VII. CONCLUSIONS 

As turbulent flow evolves, every so often we catch a 
glimpse of a familiar structure. For any finite spatial reso- 
lution and time, the flow follows unstable coherent struc- 
tures belonging to an alphabet of representative states, 
here called 'templates'. However, in the presence of sym- 
metries, near recurrences can be identified only if shifted 
both in time and space. 
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In the method of sections (along time direction) and 
slices (along spatial symmetry directions), the identifi- 
cation of physically nearby states is achieved by cutting 
group orbits with a finite set of hyperplanes, one for each 
continuous parameter, with each time trajectory and 
group orbit of symmetry-equivalent points represented 
by a single point. The method of slices is akin to (but 
distinct from) cutting across trajectories by means of sec- 
tions. Both methods reduce continuous symmetries: one 
sections the continuous-time trajectories, the other slices 
the layers of the onion formed by group-orbits. Both are 
triggered by analogous conditions: oriented piercing of 
the section and oriented piercing of the slice. Just as a 
Poincare section goes bad, the slice hyperplane goes bad 
the moment trans versality is lost. A slice, however, is 
emphatically not a Poincare section: as the first step in 
a reduction of dynamics, a slice replaces a trajectory by a 
continuous symmetry-reduced trajectory, whereas in the 
next step a Poincare section replaces a continuous time 
trajectory by a discrete sequence of points. 

The main lesson of the visual tour undertaken above is 
that if a dynamical problem has a continuous symmetry, 
the symmetry must be reduced before any detailed anal- 
ysis of the flow's state space geometry can take place. 
So far, this has only been achieved for transitionally tur- 
bulent numerical pipe flows, 11 resulting in the discovery 
of the first relative periodic orbits embedded in turbu- 
lence. In the future, it should be the first step in the 
analysis of any turbulent data, numerical or experimen- 
tal. Once symmetry reduction is achieved, all solutions of 
a turbulent flow can be plotted together: all symmetry- 
equivalent states are represented by a single point, fami- 
lies of solutions are mapped to a single solution, relative 
equilibria become equilibria, relative periodic orbits be- 
come periodic orbits, and most importantly, the analy- 
sis of the global dynamical system in terms of invariant 
solutions and their stable/unstable manifolds can now 
commence. 
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